Critical polymer— polymer phase separation in ternary solutions 
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We study polymer-polymer phase separation in a common good solvent by means of Monte 
Carlo simulations of the bond-fluctuation model. Below a critical, chain-length dependent concen- 
tration, no phase separation occurs. For higher concentrations, the critical demixing temperature 
scales nonlinearly with the total monomer concentration, with a power law relatively close to a 
renormalization-group prediction based on "blob" scaling arguments. We point out that earlier 
simulations and experiments have tested this power-law dependence at concentrations outside the 
validity regime of the scaling arguments. The critical amplitudes of the order parameter and the 
zero-angle scattering intensity also exhibit chain-length dependences that differ from the conven- 
tional predictions but are in excellent agreement with the renormalization-group results. In addition, 
we characterize the variation of the average coil shape upon phase separation. 

PACS numbers: 61.25.Hq, 64.75.-|-g, 82.35.Lr, 82.60.Lf 
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C/3 ; I. INTRODUCTION 

Understanding the incompatibility of chemically different polymer species is a topic of great practical importance 
that consequently has received widespread attention over the last decades. Although a very large number of ex- 
I periments has been performed — including determination of the coexistence curvei, the Flory-Huggins interaction 
. paramete rnTi^4 , and of critical exponents^*^ — several basic theoretical predictions^iSiSiiSiiiii^ have proven difficult to 
O ' verify, as they rely on specific assumptions, such as monodisperse and symmetric conditions, in which the solvent inter- 
^ O^ . acts equally with both polymer species. While, on the other hand, these conditions can easily be realized in computer 
simulatio ns, tho se mostly have focused on polymer demixing in binary blends, i.e., at relatively high concentrations 
(cf. Refs. Il3ll4|) . Sariban and Binderii^iiSiil have performed pioneering computational work on polymer-polymer 
^ ■ demixing under more dilute conditions, using a simple cubic lattice model. The fact that little additional compu- 
tational work appears to have been done since then has motivated us to further explore polymer demixing in the 
. presence of a solvent. 

■ Flory-Huggins (FH) theoryi^ predicts that two incompatible polymer species in a nonselective good solvent will 
, phase separate at arbitrary dilution, provided the temperature of the system is sufficiently low. Under the appropriate 
ly-j ■ conditions, the phase transition will be continuous and the corresponding critical temperature is predicted to 
. increase linearly with the total monomer concentration (j) over the entire concentration range. However, this mean- 
"Y^ field approach assumes homogeneous monomer densities for both species and ignores the chain connectivity. Below 
' the dilute-semidilute threshold 0*, this connectivity causes individual chains to be separate, swollen coils that do 
^ not interpenetrate. Thus, de Gennes^ has claimed that no phase separation occurs for concentrations (j) < cff . 
I ' Furthermore, in the semidilute regime, (f> > (f>* , each chain can be viewed as a succession of "blobs," where each blob 
contains several monomers. In a good solvent, a blob does not contain monomers of other chains due to the excluded- 
volume interactions. This reduces the number of contacts between monomers of different species and consequently 
^ lowers the critical temperature compared to the FH prediction. Specifically, the critical temperature no longer varies 
linearly with concentration, but decreases superlinearly with decreasing concentration^ii^ Employing renormalization- 
group (RG) techniques, Joanny and co-workersiS observed that demixing under good solvent conditions is driven by 
corrections to scaling. Schafer and Kappeler— calculated the corresponding spinodal by means of an RG approach 
5— i ' and found a nonlinear dependence of the critical temperature on concentration that differs from the prediction of 
5^ , Ref. O Broseta et ali^ subsequently extended the RG treatment to predict the chain-length dependences of critical 
amplitudes of the order parameter, zero-angle scattering intensity and the correlation length. These predictions have 
partially been tested in Ref. ITtI where the concentration dependence of the critical temperature in the semidilute 
regime and the chain-length dependence of the critical amplitude of the order parameter have been investigated. 

In this paper, we extend the work of Ref. ^3 by means of Monte Carlo simulations in which we employ the 
bond-fluctuation modeli2i2S rather than a simple cubic lattice model. This model is capable of a more realistic 
representation of mixtures of flexible chains, as it permits variation of the length of chain segments and a much 
smoother variation of the angle between successive segments. We investigate phase separation as a function of total 
monomer concentration and test the various theoretical predictions for the critical properties, where we also address 
an alternative scenario leading to a nonlinear variation of the transition temperature with concentration. In addition, 
our simulations explicitly cover the transition from the dilute to the semidilute regime in ternary solutions, which 
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requires (computationally demanding) calculations in the grand-canonical ensemble. Following our earlier work-^^, we 
also investigate the effect of this transition on shape properties of individual polymer coils. The shape variation of 
polymers upon phase separation is of interest because of its influence on the physical properties of ternary solutions. 
As was first recognized by Kuhn^^, the shape of a flexible polymer is ellipsoidal rather than spherical. This may 
affect, e.g., the flow properties of polymeric fluids^'^ and the polymer- induced depletion potential in colloid-polymer 
mixturesi^ Furthermore, Murat and Kremei»24 proposed a coarse-grained model for the study of phase separation 
of polymer blends, in which each coil is replaced by a "soft" ellipsoidal particle. Since individual segments are no 
longer modeled explicitly, this may permit the observation of the demixing process at longer times scales. However, 
it also requires accurate knowledge of the actual ellipsoidal shape of the individual chains. For polymer chains in 
the dilute limit, corresponding studies have indeed been performed. Sole and StockmayeriSSiSi first introduced the 
radius-of-gyration tensor. The three eigenvalues of this gyration tensor are the squared principal components of the 
radius of gyration. From Monte Carlo simulations of an ideal polymer chain on a cubic lattice, they found that 
these three components were very different, implying an asymmetric polymer shape and confirming Kuhn's original 
observation. The asphericity A was introduced to characterize the coil shapei22*22 It takes values between (sphere) 
and 1 (rod) and was calculated for ideal and self-avoiding chains.'-^-"^° Since this asphericity does not distinguish 
prolate-ellipsoid and oblate-ellipsoid shapes, another parameter S was introduced^^ and calculated analytically for 
ideal and self-avoiding chains. However, in order to simplify the calculations for A and S (both of which are defined 
below), expectation values of ratios were replaced by ratios of averages. This approximation was found to overestimate 
the asphericity of polymer chains. '^^■■^^ Numerical calculations -^-^-'^i have supported the analytical calculations. In 
recent years, it has actually become possible to observe the shape asphericity experimentallyi^^*^ These experiments, 
as well as the vast majority of the theoretical work, focus on the shape of a single polymer chain in a highly dilute 
solution. Studies of the influence of concentration are rare and have essentially shown that the asphericity of athermal 
chains diminishes only very gradually as a function of increasing concentration^2iS In a poor solvent, the reverse 
effect occurs, due to the coil-globule transition taking place in the polymer-lean phasei^ The dependence of polymer 
shape on both concentration and solvent quality has motivated our earlier study of the shape change of polymers 
upon phase separation in ternary solutions^ It was found that phase separation strongly influences the shape of the 
minority component in a given phase. Here, we expand this work by investigating the temperature dependence of 
both the asphericity A and the prolate-oblate parameter S at fixed total concentration. 



II. THEORETICAL BACKGROUND 



Our study consists of three main aspects: (i) The scaling of the critical temperature Tc with total monomer 
concentration and degree of polymerization N; (ii) the chain-length dependence of critical amplitudes; (iii) the 
shape variation of polymers upon phase separation. In this section, we provide the necessary theoretical background 
and review the corresponding predictions. 



A. Scaling of the critical demixing temperature with concentration and degree of polymerization 

We consider a monodisperse solution containing polymers of two different species, denoted A and B, with degree 
of polymerization Na — Nb — N. In addition, the solvent is of the same quality for both species, and their 
chemical potentials are equal. Thus, a so-called symmetric ternary system is realized, in which both species behave 
identically. The only distinction between monomers of type A and monomers of type B is their mutual repulsion of 
strength eab > 0. Identical monomers experience an attraction eaa — ebb < 0. Accordingly, phase separation occurs 
at sufficiently low temperatures. Owing to the symmetric properties of both species, this transition is continuous if the 
chemical potential difference vanishes. FH theory, assuming complete screening of the excluded-volume interactions, 
predicts a linear relation between the corresponding critical temperature Tc and the monomer concentration (/) (cf. 
Ref.i), 

oc , (1) 



where Ae = (2£ab — £aa — eBB)/2. Equation (^) implies that polymer-polymer (PP) phase separation occurs for 
arbitrarily low (/>, provided that the temperature is sufficiently low. However, in the dilute regime (0 < 0*), the 
excluded-volume interaction is not screened at all and polymer coils essentially do not interact. Accordingly, de 
Gennes^iiS argued that no PP phase separation can occur for concentrations below 0*. At the overlap threshold (p* 
phase separation can be induced, provided that the segregation factor is sufficiently strong (or the temperature 
sufficiently low). In the semidilute regime (0 > (p*), the tendency to phase separate is suppressed as well, and 
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demixing will take place at a lower temperature than predicted by FH theory, due to the partial screening of the 
excluded-volume interaction. According to blob scaling arguments^, in this regime each chain can be viewed as a 
succession of blobs. Each blob only contains monomers of a single species, so that the number of A-B contacts is 
reduced and phase separation is correspondingly suppressed. Indeed, de Gennes has predicted the critical temperature 
to scale as2ii& 

^ odr ^ =^ -r'-'' i<i>>r)- (2) 

Alternatively, phase separation can be induced by varying the concentration at a fixed temperature T. At low 
temperatures, it follows from the previous argument that the total monomer concentration must be increased to (jj* . 
However, if the segregation factor is weak (i.e., at high temperatures), the homogeneous phase persists into the 
semidilute regime and separation only occurs at a critical concentration 

(/-c oc (r/7V)3--i ~ (r/7V)0-764 (0>0*), (3) 

where v = 0.588 is the scaling exponent for the end-to-end distance R (or the radius of gyration Rg) as a function of 
N in dilute solution. We observe that 4>c exhibits the same chain-length dependence as the overlap threshold (j)* . By 
contrast, an RG approach^*^ indicated that for the fixed point corresponding to phase separation under good solvent 
conditions, the chemical mismatch between A and B is an irrelevant parameter. Indeed, phase separation is predicted 
to be driven by the corrections to scaling. Consequently, the scaling predictions Eqs. (|2J) and ||2J) are modified, as 
shown by an RG calculation of the spinodal^ Specifically, it was found that for demixing in the semidilute regime, 
at a temperature T, the critical concentration scales as, 

0c oc (T/7V)(3'^-i)/(i+^) ~ (r/7V)"-624 . (4) 

From the analogy between a polymer solution and the n- vector model in the limit n — > 0, the exponent x can be related 
to the negative of the crossover exponent of the isotropic fixed point of the n-vector model with cubic anisotropy^SiiS 
This crossover exponent has been calculated in a third-order e-expansion^'', and from a Pade-Borel resummation 
the value x — 0.225 (5) was obtained.^^' The most notable difference between Eqs. Q and ^ is that in the RG 
result the critical concentration (at fixed temperature) decreases with chain length at a slower rate than the overlap 
threshold 0* ~ l/iV'^'^"^, i.e., for sufficiently long chains (and weak segregation factor) phase separation sets in at 
higher concentrations than predicted by Eq. For experiments at fixed concentration [cf. Eqs. and the 
critical temperature is now given by 

Interestingly, Olvera de la Cruz^^ generalized the work of Ref. to microphase separation in diblock copoly- 
mer solutions and predicted that the order-disorder transition temperature Tqdt exhibits the same concentration 
dependence as described in Eq. |(SJ). Lodge et al.^^ found very good agreement with this prediction for solutions of 
poly(styrene-6-isoprene) (PS-PI) diblock copolymers, whereas the transition temperatures for poly(ethylenepropylene- 
&-ethylethylene) solutions followed a different power law. While the results for PS-PI solutions can be viewed as a 
confirmation of the RG scenario, the agreement between theory and experiment may well be fortuitous, as most 
data were taken in the concentrated regime where the blob scaling approach is invalid. Guenza and Schweizes^ 
subsequently studied the order-disorder transition by means of PRISM theory and observed that local concentration 
fiuctuation effects also imply a nonlinear variation of Todt, resulting from a concentration-dependent local correlation 
hole affecting the x parameter. 



1 1 1 

(X — 4> ^I'-i ~ — f 



(6) 



Thus, while the integral equation theory cannot quantitatively capture the critical fiuctuations that lead to the RG 
result, it predicts a numerically similar power law that results from a very different mechanism. Unlike the blob 
scaling arguments, this mechanism explicitly applies to the concentrated solution regime. 



B. Scaling of critical amplitudes with degree of polymerization 



Critical properties only exhibit nonclassical scaling behavior within a certain region around the critical point, 
determined by the Ginzburg criterion!^ For larger deviations from Tc, i.e., when \t\ 3> G, where t is defined as 
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(T — Tr\/T r and G denotes the Ginzburg number, classical or mean-field-like critical exponents are observed (cf. 
Rcfs. lisEo) . For polymer mixtures, the Ginzburg number decreases as 1/-/V, so that Ising-like critical exponents 
are only observed in a very narrow temperature region around the critical point According to Broseta et ali^^, 
the modified chain-length dependence of the critical temperature and critical concentration described in Sec III Al 
also affects the Ginzburg criterion, so that the number of blobs per chain N cx (^c/fl!)*)^/^^''"^) - replaces 
N . Thus, nonclassical behavior is observed within the considerably larger region |t|7V^/(i+^) < 1. In addition, the 
chain-length dependence of all critical amplitudes is modified. In this work, we test these predictions for the order 
parameter as well as for the zero-angle scattering intensity, for which — to our knowledge — it has not been verified 
before. In our (semi-)grand-canonical simulations (see Sec. IIII|) . the order parameter m is defined as 

- riB , , 

TO = ■ , (7) 

where ua, tib are the numbers of A and B polymers, respectively. This corresponds to the concentration difference (for 
either species A or species B) between the A-rich and the B-rich phases. In simulations of finite systems, the ensemble 
average (m) vanishes^, which is resolved by employing the absolute value of the order parameter, (|to|). For the sake 
of consistency, we thus phrase the following theoretical expressions in terms of (|to|). Composition fluctuations are 
probed via the zero-angle scattering intensity (osmotic compressibility)'*^ 

5cou(0) = nA^0((|m|2) - {\m\)^) oc (b^V{{\m\^) - i\m\f) , (8) 

where n = riA + riB is the total number of polymers. Near criticality, both (|to|) and 5*0011(0) exhibit a power-law 
dependence on the reduced temperature i, 

(|to|) = Bmtf , (9) 

5coii(0) = f(iV)|i|-'^ , (10) 

where (3 and 7 denote critical exponents and B{N) and r(A^) are A^-dependent amplitudes. Within the asymptotic 
scaling region, f3 and 7 assume their Ising values, (3 ~ 0.327 and 7 ~ 1.237^ The chain-length dependence of the 
amplitudes follows directly from the observation that t/G — tN constitutes the proper scaling variable and that B 
and r must be independent of N in the mean-field limit. Thus, 

B{N) oc N^'-^/^ ~iV-0-i73^ (11) 
f(A^) cx N ■ N^-'^ ~ iV°-'^63 ^ (j^2) 

where r(A^) involves an additional factor N because of the prefactor in Eq. (|SJ. The revised scaling of B{N) and 
f (TV) with N yieldsiSiil 

B{N) oc ^ ^x(/J-l/2)/(l+x) ^ ^-0.0318 ^ (;^3) 

f (A^) cx N ■ N^-'' = ivi+^(i-7)/(i+:';) ^ Tv^-^s^ . (14) 

Equation H13|l was verified in Ref. ^3 It is one of the goals of the present work to test the prediction Eq. H14|l for the 
compressibility and to reproduce the scaling for the order parameter in the context of the bond fluctuation model. 
A flnal noteworthy point concerns the so-called Fisher renormalization of critical exponents^, which is expected to 
occur for experiments and simulations at fixed concentration. It has been predictedi^ that for polymer demixing 
this renormalization occurs within an exceedingly narrow range around the critical temperature, whereas regular 
Ising-type exponents are predicted to occur outside this range (but within the nonclassical regime predicted by the 
Ginzburg criterion). 

C. Shape variation of polymer coils upon phase separation 

The ellipsoidal shape of a polymer coil is characterized by the eigenvalues Ai ^ A2 ^ A3 of the radius-of-gyration 
tensor Q, which is defined asSLffi 

1 ^ 

'^"/^ = ^772 H t''^" ~ '"j,"]K./3 - ''i,/?] ' (15) 
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where represents the position of the ith monomer along the chain and a, P = 1,2,3 denote cartesian components. 
An important measure is the asphericity 

which takes values between (sphere) and 1 (rod). In dilute solution it approaches a universal value as N —^ oo, 
estimated as 0.415 from first-order e-expansions and as 0.431 from simulationsi^Si^iS In the melt limit, where the 
chains behave ideally, this value is anticipated to decrease to the (exactly known) RW value 0.39427. . ..^^ Another 
measure of the polymer shape is 5*, 

s ^ (17) 

where A is the average value of Ai, A2 and A3. S distinguishes oblate-shaped polymers {—^<S<0) from prolate- 
shaped chains (0 < 5 < 2). In the dilute limit it approaches a universal value 0.541, as estimated from simulations 
Upon approach of the melt limit it decreases to a value that is estimated as 0.475 from numerical integration^^ and 
0.478 from simulationsi^^ 



III. SIMULATION MODEL AND DETAILS 



The bond fluctuation model (BFM) has been introduced for the Monte Carlo simulation of polymer systems in 
Refs. and In this lattice model, a chain consists of connected units, each of which occupies a cubic lattice 
cell. Every unit represents a Kuhn segment corresponding to 3 to 5 real monomers^ The units cannot overlap and 
are connected by bond vectors with lengths between 2 and \/T0 lattice constants; restricting the vectors to this set 
prevents the crossing of bonds. Important differences between the BFM and a self-avoiding random-walk (Verdier- 
Stockmayer— ) model are the flexibility of the segment length (the distance between two connected units) and the 
resulting increase in the number of possible bond angles (i.e., the angle between two adjacent segments). Monomer 
interactions are implemented by means of a square-well potential with a range of \/6 lattice constants, which covers 
54 out of the 108 possible positions for neighboring segments. The attractive coupling strength between segments 
of the same species is set to eaa — £bb = — l/fcH?"- Unlike species repel each other with a strength eab — 
where (5 > is a variable parameter. Throughout this paper, units are chosen such that — 1 and temperatures 
correspond to the inverse couphng constant S/eab — — I/^aa- The solvent is represented by empty lattice sites, and 
the polymer-solvent interaction eas and the solvent-solvent interaction ess both vanish. 

The ternary system described here can undergo both polymer-polymer (PP) and polymer-solvent (PS) demix- 
ing. Although the former typically occurs at a higher temperature than the latter, it follows from the predictions 
of Sec. ^ that the PP demixing temperature decreases with decreasing concentration. Indeed, we found that for 
simulations using the above-mentioned interaction parameters PP demixing is preempted by PS separation near the 
dilute-semidilute threshold <j)* . In order to be able to determine the concentration dependence of the PP demixing 
temperature, we therefore explicitly suppress PS phase separation in a subset of our calculations. The data presented 
in Sections IIV Bl and IIV Dl fwith the exception of Fig. I14|) pertain to eaa — £bb = 0. This eliminates all attractive 
interactions and, apart from the excluded- volume interactions, only leaves the A-B repulsion nonzero. 

All simulations are performed for symmetric, monodisperse systems (TVa = = N) on simple cubic lattices 
with linear dimension L and periodic boundary conditions. In order to permit a finite-size scaling analysis for the 
determination of critical properties, four different values for L are studied. Most simulations are performed in the 
semi-grand-canonical ensemble^'"^, in which the total monomer concentration — (l)A + (f>B is kept constant and only the 
identity of chains is changed, governed by their interactions in a given configuration as well as their chemical potential 
difference. The symmetry of the system makes such Monte Carlo moves possible and also guarantees that the critical 
point will occur for identical chemical potentials. The A-B changes are supplemented by local (translational) segment 
moves and reptation-like moves. For strong repulsions, polymer demixing occurs upon variation of concentration 
rather than temperature (this is the predicted demixing near (/> = (/)*, see Sec. Ill A)l . Therefore, these simulations are 
carried out in the full grand-canonical ensemble, in which the total monomer concentration fluctuates as well. These 
calculations, which are particularly computationally intensive, are performed using a variant of the recoil-growth 
schemei^SiS The chemical potential of both species is varied, but the chemical potential difference is fixed at zero in 
order to maintain critical demixing. 

In the semi-grand-canonical simulations, all properties are sampled every 50 or 200 sweeps (depending on </> and T), 
where a sweep corresponds to a sequence of, on average, three reptation moves per chain, one local move per monomer 
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TABLE I: Summary of linear system sizes L studied for different chain lengths A*' and concentrations (j>. Each system contains 
n = <f)L^/{?,N) chains. 
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FIG. 1: Finite-size scaling behavior of the characteristic temperatures T{L) determined from crossing points of the fourth- 
order amplitude ratio Ql and from maxima of the specific heat Cv- Data apply to a ternary polymer solution with chain 
length A'' = 20, total monomer concentration 4> = 0.32 and repulsion parameter 5 = 1. Both cases extrapolate to virtually 
identical estimates for the critical demixing temperature Tc. 



and a semi-grand-canonical move for one quarter of all chains. After equilibration, which comprises 4 000 samples, 
100 000 samples are obtained for each state point. For the grand-canonical simulations, properties are sampled every 
sweep, which corresponds to 200 attempts to insert or delete a polymer chain followed by a semi-grand-canonical 
move for half of all the chains. In this case, equilibration corresponds to 500 sweeps and each production run to 
20 000 to 60 000 sweeps. The data are analyzed by means of multiple-histogram reweightingtSS For the semi-grand- 
canonical data only temperature reweighting is performed, whereas the grand-canonical data also permit rewcighting 
with respect to the total monomer concentration (cf. Figs.|3^ and^. 



IV. SIMULATION RESULTS 



A. Scaling of the critical demixing temperature in semidilute solutions 

In order to determine the critical properties, we first locate the critical temperature for the polymer-polymer 
phase separation, as a function of concentration (j) and chain length N . For N = 10, 20 and 40 we have simulated 
cubic cells containing up to 1792, 2000 and 750 chains, respectively. TablcQllists the system sizes and concentrations 
employed. We locate by means of the finite-size scaling properties of the fourth-order amplitude ratio Q 

Ql - {mY/{m^) (18) 

and the specific heat Cv{L), 

Cv{L) - {{E^) - {Ef)/kBT^ . (19) 

For each pair of system sizes Li and L2, the curves for Ql as a function of temperature exhibit a crossing point that 
defines a characteristic temperature Tq{L), where we choose L = {Li + L2)/2. Likewise, Cv{L) exhibits a maximum 
at a temperature Tc{L). Both characteristic temperatures approach the critical temperature in the thermodynamic 
limit. Figure n ihustrates this for a typical system with N = 20, (p = 0.32 and 6=1. In this figure, we also exploit 
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FIG. 2: Scaling of the inverse critical demixing temperature N/Tc with concentration (j>, plotted on a double-logarithmic scale. 
The line, which represents a fit to the data for N = 10, describes a power-law dependence Tc oc (f)'^'^'^^. The open triangles 
refer to data for a stronger repulsion S between unlike monomers and hence lie at systematically higher temperatures. 



the theoretically predicted leading finite-size scaling behavior of the characteristic temperatures, 

1 1 D 



T{L) Tc Li/i^i 



oo , (20) 



where the coefficient D is nonuniversal and also depends on the thermodynamic property for which T{L) is determined. 
Vi denotes the critical exponent for the correlation length and takes the Ising value 0.630;— Extrapolation of Tq{L) 
and Tc{L) using Eq. H20|l yields virtually identical estimates for Tc. 

Having obtained the critical temperatures for the 15 cases listed in Table HI we investigate their scaling with chain 
length and concentration. As shown in Fig. |21 for fixed chain length and monomer repulsion ^, the inverse critical 
temperatures exhibit a power-law dependence on 0, 



N 1 
T 



« - ■ (21) 



The results for A'' = 10 can be described by an exponent k = 1.808 (5), whereas least-square fits for iV = 20 and 
N = AQ lead to fc = 1.870 (3) and k = 1.880 (4), respectively. This clearly refutes the mean-field result Eq. 1^1, 
according to which Tc increases linearly with concentration {k = 1). On the other hand, we note that our result also 
shows a stronger variation with than the nonlinear relation Eq. The blob-scaling arguments underlying the 
theoretical prediction rely on the assumption that the radius of each blob scales as J^^iob' where nbiob is the number 
of monomers in a single blob. It would have been surprising if this scaling behavior would be obeyed for the rather 
short chain lengths employed here, which necessarily contain only few, relatively small blobs. In particular, for small 
blobs the excluded- volume interactions are partially screened, reducing the magnitude of the effective exponent v^s- 
Indeed, the observed values for k = 1.80-1.88 would correspond to an effective blob scahng exponent fes = 0.55-0.56, 
i.e., only a rather small deviation from u = 0.588. However, we note that our findings are also compatible with the 
alternative expression Eq. for an effective scaling exponent v^s — 0.54-0.57. 

Interestingly, Sariban et ali^, employing short chains and a simple cubic lattice model, appear to have found 
excellent agreement with a power law k = 1.60, precisely matching the blob-scaling result Eq. jSJ- Careful examination 
of Fig. 9 in Ref. 0, however, shows that all data that support this scaling result were obtained for high monomer 
concentrations, cj) = 0.40, 0.60, 0.80 and chain lengths N ^ 32. For such concentrated solutions, the blob scaling 
arguments certainly do not apply, casting a doubt on the origin of the observed power-law behavior. Inclusion of the 
data for (f> — 0.20 leads to exponents that are significantly larger in magnitude. We note that even for the same value 
of the scaling variable employed in Ref. T? [TV^-'^/^^""^^], the data for the longest chains {N — 64) lead to a power law 
that differs appreciably from k — 1.60. This brings the observations in Ref. [13 considerably closer to our data for the 
bond fluctuation model. All data certainly appear consistent with a slow approach of the predicted limiting behavior 
Eq. jSJ, but it also cannot be excluded that simultaneously an alternative mechanism^ applies that leads to a rather 
comparable power-law behavior in the concentrated solution regime. 

Another aspect in Fig. [21 is the dependence of the critical demixing temperature on the chain length. Although the 
plotted data for S = 1, in which Tc is scaled by N, do not exactly coincide, the remaining chain length dependence is 
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FIG. 3: (a) Zero-angle scattering intensity Scoii(O) (scaled by the chain length A'^) as a function of monomer concentration (j> 
for a ternary solution containing two polymer species with degree of polymerization A'^ = 20. The data are obtained from 
grand-canonical simulations at fixed (reduced) temperature T = 0.05. Away from the maximum the data for different linear 
system sizes L converge rapidly, whereas the maximum itself exhibits strong finite-size efi'ects. (b) The height of the maximum 
increases with system size as U^^"'^ ~ L^'^^, where 7 and v\ are critical exponents belonging to the Ising universality classS 
This confirms the occurrence of a critical phase transition as a function of (j). As discussed in the text, polymer-solvent 
separation is suppressed by eliminating attractive interactions between monomers of the same species, £aa = ebb = 0. 



weak and suggestive of an additive small- correction that vanishes in the limit N ^ 00. In order to show the effect 
of the A-B repulsion on Tc, we also include data for S = 2 ioi N — 20. These results exhibit the same power-law 
dependence on concentration, but phase separation occurs at systematically higher temperatures than for 5 = 1, as 
expected from the stronger A-B repulsion. 



B. Polymer— polymer phase separation and the role of the dilute— semidilute threshold 

According to FH theoryi^, PP phase separation occurs at arbitrarily low concentration. On the other hand, de 
Gennes:^'— has predicted that dilute solutions {(j) < (j>*) exhibit no phase separation, whereas a critical demixing 
transition takes place for symmetric systems if the A-B repulsion is strong enough (or temperature sufficiently low). 

In order to test this scenario, we consider a ternary solution with TVa = Nb = TV = 20 at very low, fixed temperature 
T — 0.05 (strong repulsion regime). We perform simulations in the grand-canonical ensemble at several chemical 
potentials that correspond to a range in monomer concentration and we specifically monitor the zero-angle scattering 
intensity S'coii(O), Eq. (jHJ, and the fourth-order amplitude ratio Ql, Eq. (|18|) . Since the simulations are carried out 
for four different system sizes, L = 40, 50, 60, 70, a critical transition will manifest itself via finite-size effects. The 
compressibility will exhibit a maximum that increases with system size according to a well-defined power law and 
the amplitude ratio will exhibit a universal crossing point. The compressibility data shown in Fig. which were 
obtained by means of histogram reweighting, exhibit precisely this behavior. For each system size, S'coii(O) displays a 
maximum S^if{0) as a function of monomer concentration. This maximum scales with L as L^/"^! ~ L^-^^ (Fig.lSt). 
The power law is characteristic for the compressibility at a critical phase transition with a one-component order 
parameter. We note that these simulations are performed at very low temperatures. Nevertheless, phase separation 
only occurs when a certain critical concentration has been reached. The fourth-order amplitude ratio of moments 
of the order-parameter distribution confirms these observations. As shown in Fig. ^ the curves for Ql for different 
system sizes exhibit a crossing point at a concentration that is close to the concentration of maximum compressibility 
(Fig. The crossing curves are very similar to those commonly employed to determine critical temperatures; the 
curves for successive system sizes do not all cross at a single concentration owing to finite-size corrections. However, 
it is noteworthy that the value of Ql at the crossing point approaches the universal Ising value Ql = 0.6233 (cf. 
Ref. ^Ed) with increasing L. 

Tabic HTl lists the critical concentration <j>c determined from both the compressibility and the fourth-order amplitude 
ratio for chain lengths 10, 20, and 40. As an alternative approach, we exploit the notion that, at fixed concentration (j>, 
the curves for Ql{T) for two different system sizes L must exhibit a crossing point, if the system undergoes phase 
separation. Thus, if such a crossing point is observed, 4> > 4>c and if no crossing point can be detected, (j) < 4'c- 
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FIG. 4: Fourth-order amplitude ratio Ql as a function of monomer concentration for different system sizes L. The data pertain 
to the system described in the caption of Fig. |^ The crossing points, which occur near the location of the compressibihty 
maxima in Fig.|^, approach the concentration at which critical polymer-polymer demixing takes place. 



TABLE II; Critical concentration for phase separation at low temperatures (strong segregation factor), is determined from 
extrapolation of the maxima of the compressibility and of the crossing points of the fourth-order amplitude ratio Ql (cf. Figs. 
l^and^J- 4> is determined by detecting whether, at fixed concentration, two curves for Ql exhibit a crossing point at some 
temperature. 



N 


4>c 





10 


0.131 ±0.005 


0.13375 ±0.00125 


20 


0.101 ±0.005 


0.10375 ±0.00125 


40 


0.077 ± 0.003 


0.07625 ± 0.00125 



The resulting estimates {(f) in Table HTjl indeed agree reasonably well with those for 0c- Since the second method is 
computationally more efficient , we have used it to estimate the critical concentration for chain lengths as large as 
N = 320, see Fig.El 

According to de Gennes^, in a good solvent the overlap threshold scales as 

cf)* cx N/Rl oc N^-^" = iV-°-^^4 , (22) 

and for the low-temperature regime 4>c is predicted to follow the same scaling behavior. Our findings for 0c (Fig. IS)) 
exhibit a slowly varying, effective power-law behavior, which for the longest chains approaches 0c oc TV^O-SO. ft is well 
possible that the trend in the effective exponent continues and reaches the behavior predicted in Eq. (|22|l . We note 
that the analysis of Broseta et ali^, which modifies the scaling of 0c from Eq. ^ to Eq. Q), only appears to apply to 
the semidilute regime. This would imply different scaling behavior for the weak and the strong segregation regimes, 
and it is therefore not entirely clear whether de Gennes' original prediction for the low-temperature regime (0c — 0*) 
indeed remains unaltered. 

Finally, the critical concentrations for the low-temperature regime can be combined with the results for the systems 
listed in Table ^ to assemble a phase diagram showing the critical lines for three different chain lengths in the 
concentration-temperature plane, see Fig. |H| The critical lines in the semidilute regime are terminated by the lines 
labeled 0*. These cutoff lines are drawn vertically to reflect the observation that in this model 0* is insensitive to 
temperature variation. 

C. Critical amplitudes in ternary solutions 
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Near a critical point, finite-size scaling theory implies 
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FIG. 5: Scaling of the critical concentration (j>c with A'^ on a log-log scale. All results pertain to a system with eaa = ebb ~ 
in the strong segregation regime {1/T — 20). The estimates for N > 40 were obtained from the detection of crossing curves 
for Ql{T), as explained in the text (cf. in Table [HJ. The dashed line represents the theoretical prediction that the critical 
concentration in the strong segregation regime is proportional to the overlap threshold. The solid line represents the effective 
power-law behavior observed for the longest chains investigated in this work (120 ^ A'' ^ 320). 
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FIG. 6: Phase diagram for polymer-polymer separation as a function of total monomer concentration and temperature, for 
different A''. All data pertain to systems without monomer attractions, £aa = £bb ~ 0, and hence differ from those presented in 
Fig. |5| The curves (which are drawn as guides to the eye) separate the mixing region (left-hand side) from the demixing region 
(right-hand side). The curve b indicates the weak segregation regime, in which phase separation occurs at high concentrations, 
whereas the dotted line a reflects the strong segregation regime where phase separation takes place for concentrations above 
the overlap threshold. (For clarity, only the curves for iV = 40 are labeled.) 



and 

S,on{0)L-^/'''/cp^ = S{tL^/'''), (24) 

where / and S are universal scaling functionsiSS In the finite-size scaling limit, tL^/'^^ ^ 1, these equations yield 
the finite-size scaling behavior at the critical point, commonly employed for the numerical determination of critical 
exponents^ Outside the finite-size scaling regime, tL^^'^^ ^ 1, but still sufficiently close to the critical temperature 
to be within the scaling regime t ^ 1, the scaling functions reproduce the critical behavior in the thermodynamic 
limit, i.e., /(x) ^ and S{x) ^ x^^ . Here, we focus on this second regime. Assuming Ising values for the exponents, 
we plot {\m\)L^/'^^ in Fig. [7|as a function of \t\L^/'^^ for a solution with = 40 and total concentration = 0.24. 
For this system, which has a critical temperature = 22.0 (1), two sets of data are plotted, obtained at T = 18.182 
and T = 20.000, respectively. Each set contains four different system sizes L, and all data turn out to collapse on a 
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FIG. 7: Finite-size scaling of the demixing order parameter (|m|) for a solution with <j) — 0.24 and A'^ = 40, near its critical 
temperature Tc — 22.0 (1). The two data sets pertain to temperatures T — 18.182 and T — 20.000, respectively. For each 
temperature, the data points correspond to system sizes L = 40, 60, 80, and 100. All points are described by a single power 
law {\t\L^/''^f, with f3 = 0.327. 
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FIG. 8: Finite-size scaling behavior of the zero-angle scattering intensity Scoii{0)/4>'^ for the system described in the caption 
of Fig. 13 Panel (a) employs the regular reduced temperature \{T — Tc)/Tc\, which does not permit a description of all 
data by a single power law, presumably because of corrections to scaling. In panel (b), the alternative temperature variable 
[exp(A'^/r) — exp(7V/rc)] is used, which indeed greatly improves the data collapse and permits the description of all data by a 
single power law {\t\L^^'^''-)~"' , with 7 — 1.237. The left-most data point in panel (b) deviates because of finite-size corrections. 



single line which, on this double-logarithmic scale, has a slope 0.327, the critical exponent /? in the Ising universality 
class. This confirms Eq. (|23|) and our assumption of Ising exponents. Likewise, data for the zero-angle scattering 
intensity S'coii/0^ are plotted in Fig. The data do not fall onto a single curve as well as the data for the order 
parameter do, which can possibly be ascribed to the fact that the reduced temperature for our data is relatively 
large and hence corrections to scaling start to become important. Therefore, we replot the same data in Fig.|S|D as a 
function of an alternative temperature variable [exp(iV/T) — exp(A^/Tc)], proposed in Ref. fisl To leading order, this 
variable equals — [(A^/Tc) exp(7V/Tc)]i. Note that, as implied by Eq. |(2U, the prefactor is independent of N. This 
variable indeed improves the scaling behavior and the data in Fig. |SJd are well described by a power law {tL^/"^)-'^, 
with 7 — 1.237. Similar scaling analyses have been carried for our results for TV = 10 and iV = 20 (not shown). 

By combining the data for different chain lengths, we can test the modified scaling of B{N) and T{N) with N as 
proposed by Broseta, i.e., Eqs. H13|l and H14|) . Figure El constitutes the counterpart of Fig. [7| displaying the scaling 
of the order parameter (|m|) outside the finite-size scaling regime (but within the critical region). In Fig. the 
conventional scaling with N is adopted, in which the reduced temperature t is replaced by tN cx t/G and (|to|) is 
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FIG. 9: Scaling of the near-critical order parameter (|m|) with system size L and degree of polymerization A'^. Panel (a) 
demonstrates that the conventional scaling [Eq. (llH l is not applicable in the semidilute regime. In panel (b), a modified scaling 
of the critical amplitude is adopted [Eq. 1131 1. which leads to an excellent collapse of all data points onto a single power law with 
Ising exponent (3 — 0.327. Data points pertain to the following systems (all at concentration tj) = 0.24): (i) A'^ = 10 (Tc = 6.309): 
T = 5.952 and T = 5.650. For each temperature, system sizes L = 30, 40, 50, and 60 are employed, (ii) iV = 20 (Tc = 11.547): 
T = 10.417 and T = 9.434. For each temperature, system sizes L = 40, 60, 80, and 100 are employed, (iii) TV = 40 (Tc = 22.0): 
T = 20.000 and T = 18.182. For each temperature, system sizes L = 40, 60, 80, and 100 are employed. 



multiplied by N^l"^ , following Eq. (|ll|l . Evidently, this does not properly describe the scaling of the critical amplitude 
with the degree of polymerization. However, replacing N by iV^/(^+^\ i.e., plotting (|m|)iV^/[^(-'^+^^lL^/'^i as a function 
of \t\N^ I '^^^^^ L}-I (Fig.Eb) leads to an excellent collapse of all data points onto a single line that describes a power 
law with exponent /? = 0.327. Thus, the order parameter scales as 



which can be simplified to 



(|m|) cx iV-(^-"-5)/(i+-)t/3 , 



(25) 



(26) 



in agreement with Eq. (|13|) . 

In a similar fashion we test Eq. (|14|l . Based upon Fig. |Slwe use the modified reduced temperature [exp(A^/r) — 
exp(A^/rc)] instead of — i, which does not affect the scaling with N . Since all data points apply to the same concen- 
tration, we investigate the scaling of 5*0011(0) / {(pN). Figure ITUl demonstrates that all data points except for those with 
the smallest values of \t\L^^'^^ (which cross over to a horizontal line representing the critical finite-size amplitude) are 
described by a power law with exponent —7 = —1.237. Thus, the following scaling behavior is recovered. 



'5'coll(0) j^-y/uiJ^T-x/(l + x) 
(f)N 



which can be simplified to 



<f)N 



cx 



^a:(l-7)/(l-|-a;)^-7 



(27) 



(28) 



consistent with Eq. (|14|l . 

Thus, we conclude that the chain- length dependence of the critical amplitudes in the semidilute regime can not 
be described by the conventional scaling laws, but that the modified scaling behavior proposed by Broseta et alm^ 
provides an excellent description already for relatively short chains. 



D. Shape variation of polymers upon polymer—polymer phase separation 

In our earlier communication^i, we studied the shape change of polymers upon phase separation in terms of 
the asphericity A (see Eq. (|16ll ). We confined ourselves to isothermal variation of the concentration in the strong 
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FIG. 10: Scaling of the zero-angle scattering intensity near criticality. Data are assembled for all systems listed in the caption 
of Fig. El Except for some data points that cross over into the finite-size scaling regime (left-hand side of the plot), all data 
are described by a single power law with exponent —1.237. This confirms the modified scaling with degree of polymerization, 
see Eq. (QIJ . 



segregation regime, i.e., phase separation near (p = (p* . This corresponds to an isotherm that intersects the line a 
in Fig. Along such an isotherm, the asphericity of the majority component in a given phase decreases merely 
slightly at sufficiently high concentration due to the screening of the excluded- volume interactions (thus, the shape 
of a coil becomes slightly more spherical upon increase of the concentration). By contrast, the minority component 
in each phase exhibits a strong variation of the asphericity. Since a typical coil belonging to the minority species is 
surrounded by polymers of the majority species, the strong repulsion "squeezes" the minority coil to a much more 
spherical shape, as confirmed by a rapid drop in A for the minority polymers. Here, we extend this work by studying 
other regions of the phase diagram, namely (i) variation of the total monomer concentration at higher temperatures 
(such that phase separation takes place in the semidilute regime) and (ii) temperature variation at fixed concentration 
(for (j) > (j)*), where we can also verify any influence of attractions between monomers of the same species (in R.ef. I21I 
such interactions were explicitly ruled out, in order to prevent polymer-solvent phase separation near the overlap 
concentration). In addition, we extend our analysis of coil shapes by means of the parameter S [Eq. (|17|) ]. which 
distinguishes between prolate-ellipsoidal and oblate-ellipsoidal shapes. 

Figure [TTk shows the variation of A for a symmetric mixture with = 20 along the T — 2.14 isotherm, which 
crosses the critical line in the semidilute regime (cf. Fig. [SJ. Just as observed for the strong segregation regime—, 
the asphericity for the majority component decreases slowly while for the minority component it drops rapidly. The 
location of the bifurcation in A is in good agreement with the corresponding critical concentration in Fig. |B| The 
nature of the aspherical shape is characterized further by the variation of the oblate-prolate parameter S over this 
concentration range (Fig. Illb ). As the concentration increases in the dilute regime, the coil shape becomes less 
prolate, until the concentration has reached its critical value and phase separation occurs. Beyond (j)c, the minority 
component (e.g., a chain of type B in the A- rich phase) then becomes less prolate at an even higher rate. On the 
other hand, the majority component becomes less spherical and more prolate immediately after phase separation. 
As explained in Ref. l2lL we ascribe this to the diminished repulsion that a typical majority chain experiences in a 
homogeneous (unmixed) phase. Although the concentration dependence of A and S is very similar, there are significant 
differences in the probability distributions for these quantities. While the distributions for A are relatively broad—, 
the distributions for S, both for the majority component (Fig. 112b .) and for the minority component (Fig. 112b ) exhibit 
a sharp peak at zero (spherical shape). For the majority component the distribution, including the peak, exhibits only 
minor variation with concentration in the dilute regime and remains virtually unchanged for (f> > (f>c- However, for the 
minority component the most significant changes in the distribution, including a significant sharpening of the peak, 
occur for concentrations in the unmixed regime (f> > 0.16. For both components, the distributions provide insight into 
the relative occurrence of various coil shapes, which can also be oblate-ellipsoidal (5 < 0). 

As evident from the phase diagram (Fig. (H)), phase separation can also be induced by lowering the temperature 
at fixed concentration. The resulting variation in coil shape is illustrated in Fig. 1131 for a system with N = 20 and 
(j) — 0.16. A bifurcation similar to Fig.lllb is observed, at a temperature that agrees with the critical temperature along 
the corresponding isochore in Fig. |S1 Unlike the behavior upon variation of the concentration, the asphericity of the 
majority component now remains constant in the unmixed phase. Since the concentration is constant, the screening of 
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FIG. 11: Characterization of the coil shape in a ternary solution, as a function of monomer concentration (j>. The temperature 
is fixed at T = 2.14 and the results pertain to a symmetric mixture with chain length A'^ — 20. Phase separation occurs 
near (j> = 0.16 and is driven purely by repulsions between unlike monomers (eaa ~ ebb ~ 0). (a) The asphericity A [Eq. lITCl l 
decreases with increasing concentration, in particular for the minority component, indicating that the coils take an increasingly 
spherical shape, (b) The oblate-prolate parameter 5* [Eq. Il7ll 1 demonstrates that the coils have a prolate shape over the entire 
concentration range. The "elongated" shape is most pronounced in the dilute regime. 
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FIG. 12: Probability distribution of the oblate-prolate parameter S of (a) the majority component and (b) the minority 
component for different concentrations (f), eit T — 2.14 for N = 20. All data pertain to eaa = £bb = 0. While the majority 
component exhibits a virtually invariant distribution for concentrations in the demixed regime i\> > 0.16, the minority component 
shows a significant sharpening of the peak at S = 0, which corresponds to a spherical coil shape. 



the excluded-volume interactions remains unchanged, and the chains only become somewhat more aspherical as phase 
separation continues, owing to the diminishing repulsion from the minority component. Since all attractions between 
monomers of the same species have been set to zero, the majority component finds itself, once phase separation is 
complete, essentially in an athermal one-component solution. 

The effect of attractions between monomers of the same species is addressed in Fig. I14L for a system with eaa = 
£bb = — 1/r. The overall behavior of the asphericity is very similar to that displayed in Fig. Illb. with a small, 
systematic lowering of the asphericity in the dilute regime, induced by the mutual attractions. Thus, the results 
presented above fFigs. [Tnand ll2fl and in Ref. |21| are not qualitatively affected by the absence of monomer-monomer 
attractions. 
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FIG. 13: Variation of the asphericity as a function of inverse temperature for fixed concentration (p — 0.16. The results pertain 
to a system with A'^ = 20 and £aa = ebb = 0. Once phase separation is complete (i.e., at low temperature) the majority 
component is essentially in an athermal one-component solution. 
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FIG. 14: Asphericity A as a function of concentration for systems with an explicit, temperature-dependent attraction between 
monomers of the same species, eaa = ebb ~ —l/T. Chain length is set to A?^ = 20 and temperature to T = 5.39. This confirms 
that the absence of A-A and B-B attractions does not qualitatively affect the results presented in Figs. Illj and ll2l 



V. CONCLUSIONS 



We have studied polymer-polymer phase separation in a common (nonselective) good solvent by means of Monte 
Carlo simulations of the bond fluctuation model. Our calculations cover both dilute and semidilute solutions, making it 
possible to distinguish critical demixing in the weak and strong segregation regimes. In the weak segregation regime 
we have determined the nonlinear concentration dependence of the critical temperature by means of semi-grand- 
canonical simulations and in the strong segregation regime, where phase separation occurs upon variation of total 
monomer concentration rather than temperature, by means of grand-canonical simulations. We observed a sudden 
drop in critical temperature near the overlap threshold, as first predicted by de Gennesi^, although the chain-length 
dependence of the corresponding critical concentration differed rather strongly from that expected for the overlap 
concentration. Also the nonlinear relation between critical temperature and critical concentration in the semidilute 
regime exhibited a power law that differs from the theoretical predictioniiiS, which we ascribe to finite chain-length 
effects. However, it is also possible that, for more concentrated solutions, similar nonlinear behavior results from a 
concentration dependence of the local correlation hole.^'^ This may aid in explaining earlier numerical results^^ and 
experiments on diblock copolymer solutions^^, although there are also quantitative differences between those results 
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and the prediction of Ref. The modified chain-length dependence of all critical amplitudes, first predicted by 
Broseta et ali^ using a renormalization-group approach, has been verified explicitly for the demixing order parameter 
and, for the first time, for the zero-angle scattering intensity (osmotic compressibility), whereas a scenario in which the 
critical amplitudes retain their unmodified chain- length dependence can be convincingly ruled out. The observation 
of unrenormalized critical exponents is consistent with the prediction that Fisher renormalization of those exponents 
only takes place within a very narrow temperature range around the critical temperature. Conversely, the observation 
of nonclassical critical exponents within a rather large temperature range is consistent with a modified Ginzburg 
criterion^^, which implies a slow crossover to classical critical exponents. 

In a preliminary report^!, we observed that phase separation causes polymer coils belonging to the minority com- 
ponent to become more spherical, due to repulsion from the surrounding polymers of the opposite species. Here, we 
have recovered this behavior for more general monomer-monomer interactions and for simulations in which phase 
separation occurs in the weak segregation regime, as well as for the situation in which phase separation is induced 
by means of temperature variation rather than variation of the total monomer concentration. In addition, we have 
characterized the shape variation more precisely by means of the distribution of the prolate-oblate parameter. 
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